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0.1 Preface 

This Thesis is a review of the inclusive decays of heavy hadron. It is based 
on the author's study of two issues which are important in the description 
of the inclusive beauty decays given by the heavy quark expansion. 

• Evaluation of the expectation values of the four-quark operators in the 
heavy quark expansion, 

• Study on the validity of the quark-hadron duality in the heavy quark 
expansion. 

We evaluate the four-quark operators in two approaches: (1) in terms 
of the oscillator strength of the hadrons in the harmonic oscillator model, 
and (2) from the difference in the total decay rates of beauty hadrons under 
SU(3) flavour symmetry. The results obtained indicate that the lifetime 
difference between B meson and A& can be explained within heavy quark 
expansion by including the third order contribution. 

To discuss the quark-hadron duality, the power corrections to the par- 
tonic decays rate are obtained in terms of the gluon mass which is assumed to 
imitate the short distance nonperturbative effects, by the choice A 2 3> A.q CD . 

The Thesis is organised in six Chapters, as follows: 

• 1 Introduction: The Standard Model is briefly discussed and then 
the Spin- flavour symmetry exhibited by heavy hadrons is reviewed. 

• 2 Four-quark operators: The author's works on the evaluation of 
four-quark operators are presented. 

• 3 Spectator effects: Presented are the calculation of ratio of life- 
times of beauty hadrons, particularly, r(A;,)/r(-B) and the inclusive 
charmless semileptonic decay rate and determination of the CKM ma- 
trix element \V u b\. 

• 4 Short distance nonperturbative corrections: The renormalons 
contribution to the heavy quark expansion is discussed based on the 
authors work. 



m 



5 Quark-hadron duality: The validity of the assumption of quark- 
hadron duality is qualitatively discussed. 

6 Conclusion: Conclusion of the Thesis with outlook of the works 
done are given. 

I am glad to acknowledge that my acquaintance with heavy quark 
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Chapter 1 

Introduction 



This Chapter consists of two sections: The first is devoted to the Standard 
Model and in the second the heavy quark symmetry and the heavy quark 
expansion are discussed n. 

1.1 Standard model 

In this section, we briefly describe the texture of the Standard Model and 
remark on the perspectives of the weak decays of heavy hadrons in the 
Standard Model. 

1.1.1 Quarks, Leptons, and Hadrons 

Today, our knowledge is, as embodied in the Standard Model, firm about 
the primitive constituents of matter pj. They are six quarks and six lep- 
tons, all fermions. The quarks are interacting through all the known three 
interactions of the elementary particles in nature: electromagnetic, weak 
and strong, but the leptons do not strongly interact. Both of them are not 
influenced by the gravity effects. The interactions are mediated by (gauge) 
bosons: photon, W ,Z° and gluon. The properties of these fermions and 



bosons are given in Tables 1.1 and 1.2 



The fermions make up the hadrons: meson, a bosonic hadron, of a quark- 
antiquark pair and baryon, a fermionic hadron, of three quarks. Quark 
structure of hadrons is described, for three quarks u, d and s, under SU(3) 
x Natural units, namely, % — c — 1, are used throughout. 



Table 1.1: Quarks and Leptons ( a in GeV/c ; in e units; c quantum num- 
bers; d charge of leptons except neutrinos is —1) 



Quark 


Mass a 


Charge 6 


Qu. nos. c 


Lepton d 


Mass 


up, u 


(1.5 - 5) .1(T 3 


+ ? 

_! 


Iz = +i 


e~ 


0.51 .10~ 3 


down, d 


(3 - 9) .1(T 3 


*.=-§ 


ft~ 


0.105 


strange, s 


0.060 - 0.170 


s = -l 


T~ 


1.777 


charm, c 


1.1 - 1.4 


i 


c = 1 


Ve 


10-15 .10" 6 


bottom, b 


4.1 - 4.4 


6 = -l 


v» 


< 0.17.10" 3 


top, t 


170 


+1 


t = 1 


V T 


< 18.10- 3 



Table 1.2: Force carriers 



Gauge Boson 


Mass 


Charge 


Force 


Strength 


photon, 7 

W 

Z 

gluon, g 




80 
90 





± 1 






electromagnetic 
weak 
weak 
strong 


1/137 
1.1663692 GeV" 2 

0.12 at M z 



flavour symmetry by the Gell-Mann - Nishijima relation: 



Q = I z + (B + S)/2 



(1.1) 



where Q is the electric charge of the hadron, in units of e, I z the third 
component of the isospin, B the baryon number, and S the strangeness 
quantum number. Each quark is assigned a colour quantum number. Each 
quark carries three colours but the quarks make hadrons colour singlet. For 
mesons, B = 0. This can be extended to the charm and beauty quarks 
as well. The situation is different in the case of top which decays before 
hadronising. The rate of top quark weak decay, T(t — ► bW) ||, is found to 
be larger than the QCD scale that corresponds to the typical hadron size. 

In the next section we present a brief review of the Standard Model 
(SM). The SM is a renormalisable field theory of the local gauge group 
SU(3)c x SU(2) L x U(1) Y §. It is a unified theory of strong (SU(3) C ), 
weak (SU(2)l) and electromagnetic (U(1)y) interactions at TeV scale. The 
SM contains two sectors: quantum chromo dynamics and electroweak H. 



1.1.2 Quantum Chromodynamics 
A prelude 

In the proposal of quark structure of hadrons made by Gell-Mann and Zweig, 
the quarks are described in terms of flavour and colour as required to explain 
the known properties of hadrons. The idea of quarks carrying colour got con- 
firmed by the measured normalisation of the decay rate r(-7r° — ► 27) and the 
cross-section o~(e + e~ — > hadrons). Using current algebra, the strong inter- 
action between quarks is shown to be mediated by massless vector bosons, 
gluons, analogous to photon of electromagnetic interaction. According to 
the colour assignment to quarks, gluons form an octet. 

Lepton (e~ jv) - proton deep inelastic scattering experiments confirmed 
the composite structure of a proton. Further, from the deep inelastic scat- 
tering experiments, it was found that quarks and gluons are freely roaming 
around inside the proton. This behaviour is explained by Bjorken through 
his scaling law that at small distance (this corresponds to large momentum 
transferred), the quarks are behaving as if they are free, but at large dis- 
tance (i. e., small momentum transferred), they are confined in pair. The 
former property is known as asymptotic freedom and the latter confinement. 

The Lagrangian of QCD 

Theory of quarks and gluons is based on the local non-abelian gauge group 
SU(3)c symmetry Q. With the work of 't Hooft on renormalisability of 
Yang- Mills theories |Q] , the theory of strong interaction is formulated. The 
basic Lagrangian is 

Cqcd = —F^FJT + J2<i La (x)^D fl q(x)l -J2 m f^ a ( x )l( x )l (1-2) 

/ / 

where the indices a, b stand for colour and / for flavour of quarks and 

F; u = d,A a u -d u A--g s f abc AlAl, (1.3) 

Dl = d; + i 9s t a Al (1.4) 

where g s is the strong coupling constant, f abc the structure constant of 
SU(3) Lie algebra and t a = ±A a , with A a being the Gell-Mann colour (3 x 3) 
matrices. In Eq. fll.2j), the first term is the gluon part, the second the 



quarks part and the quark-gluon interaction part and the third the quark 
mass contribution. The Lagrangian is invariant under the local gauge group 
SU(3). The gauge transformation is given by U(x) = e tta £ a ( x ) where £ a (x) 
is an infinitesimal gauge parameter with a running from 1 to 8. The t a are 
the generators of the gauge group. Because of the non-Abelian nature of the 
boson field, gluons can be described as self-interacting, vector bosons. The 
Lagrangian, in ( |1.2|) , is supplemented by additional terms corresponding to 
gauge- fixing, ghost and counter terms for perturbative treatment. With 
these terms, one can get the Feynman rules for QCD. 

The QCD Lagrangian, in the limit of m q = 0, where q = u,d, s, is 
invariant under the chiral symmetry group SU(3)l x SU(3)r ||. On the 
other extreme, for quarks of mass tuq 3> Aqqd, in the limit tuq — > oo, 
the Lagrangian exhibits spin-flavour symmetry Q which will be discussed 
in details in the next chapter. Besides , Cqcd possesses stability against 
weak radiative corrections due to the flavour-neutrality of the gluon. 

Perturbative QCD: Running Coupling and Running Mass 

Different renormalisation schemes are used to make the Lagrangian free of 
divergences. Then the renormalised quantities are defined at some chosen 



scale [1C]. Different choice of the scale leads to different physical values 
for the same quantity. In order to circumvent the scale dependence, renor- 
malisation group equations are used. The scale transformation, /j, — ► //, 
constitutes a renormalisation group. Under this group, the physical quan- 
tities, explicitly or implicitly scale dependent, are invariant. Such a study 



is due to Callan-Symansik. In the MS or MS scheme, for the coupling 
constant, the Callan-Symasik or renormalisation group equation is given by 

dg(n)/dlog^ = g(n)p(g(p)) (1.5) 

An expansion of (5 in g(fi) yields, at the lowest order, 

n—ikiw, A = ^-*a?> (i.6) 

M MU 16tt 2 3 3 v ' 

where Ca = 3, Tp = 1/2 and rif is the number of flavours. Now, defining 
the running coupling constant a s = g 2 /Air, with \i = Q 2 /v 2 , 

a{Q ] = A,iog(Q7A 2 )' (L7) 



A 2 = .V 41 '^) (1.8) 

where Aqcd is known as the QCD scale parameter, which is the constant 
of integration. From (|l.8| ), it is obvious that as Q 2 — > oo, a s — > 0. This be- 
haviour is the celebrated asymptotic freedom 111]] . For finite but sufficiently 
large Q 2 , a s is employable as a perturbation parameter. 
Similarly, the running mass m{Q 2 ) is given by 



m(Q 2 ) = m 




-d n 



(1.9) 



where tuq is the invariant mass and d m is known as the anomalous mass 
dimension, given by 4//3o- 

A digression 

The perturbation series of a physical quantity in terms of a s is known to few 
orders. A quantity is said to be defined completely, if the expansion in a s is 
made to all orders. But it is impossible to define them to n orders, where n is 
large. At large order, as first discovered in QED by Dyson in 1952, the series 
in a s factorially diverges. This issue has attracted interests recently after 
it was first pointed out by 't Hooft in 1977 that the inherent ambiguities in 
the perturbative definition of a quantity is O (Aqcd/Q) and beyond. These 
ambiguities are known as renormalons [p.2[ . These corrections are seen to 
be important as they would represent the first nonperturbative corrections 
to the leading order of the OPE. 

The QCD needs to be dealt with at small and large distance scales sepa- 
rately. Any physical amplitude of a QCD process described by the product 
of operators can be separated into two parts by the Wilson short distance or 
operator product expansion |plj| . As a result, the short distance physics is 
contained in the expansion coefficients, which happen to be c-numbers, and 
the local operators describe the long distance aspect. While the perturbative 
QCD successfully describes the short distance processes, the long distance 
aspects being nonperturbative are poorly understood. However, there are 



semi analytic methods such as QCD or SVZ sum rules [14] and lattice gauge 
theory [15]. 



1.1.3 Electroweak Theory 
Fermi theory of weak interaction 

With the explanation of the the radioactive decay of a nucleus, the weak 
force was discovered. For example, weak interaction drives a neutron to 
change into a proton with a lepton pair, electron and its antineutrino, inside 
the nucleus. At the quark level, it is the transition of the d quark into the 
u quark plus a lepton pair. A theoretical description of weak interaction 
emerged due to Yang and Lee, and Wu et al |l6| following Fermi. The 
Lagrangian of weak interaction is then given by 

C 

Cperrm = ~J= J"{x) jj(x) + h.C, (1.10) 

where Gf = 1.16637 ± O.OOOOlGeF -2 is the Fermi coupling constant and 
the weak current J^ is given by 

J!*(x) = P(x)^(g v - g A ~/5)n(x)e- (x) 7/ ,(l - j 5 )v e -(x), (1.11) 

where the vector coupling gy is little less than one and is given in terms of 
the Cabibbo angle to be 0.97. The ratio of the coupling constants gAJgv = 
-1.2573 ±0.0028 §. 

Towards formulating a complete theory of weak interaction, it is postu- 
lated that the weak interaction is mediated by a vector boson, massive due 
to weak interaction being point-like. The additional requirements of such a 
boson are that it has to be charged both positively and negatively because 
of charge-changing weak current and having indefinite parity to preserve the 
V — A structure of the weak current. Thus, the Lagrangian is 

Cw = -gwJ ll {x)W ll {x) + h.c. (1.12) 

where W(x) is the vector boson field and g\y the dimensionless weak coupling 
constant which is related to the Fermi coupling constant by 

g w _ G f nil) 

The Lagrangian in ( |1.12| ) is useful at low energy but not renormalisable. 
The existence of W is another way of looking at Gf. The step towards the 
renormalisable theory of weak interaction results in the unification of the 
weak and the electromagnetic interactions, electroweak interaction. 

6 



An interlude 



The current-current interaction described by ( 1.12J ) is applicable only for 
charged current. There is also a neutral current mediated weak interaction. 
If both are combined, then three massive vector bosons correspond to the 
isospin triplet of weak currents. This obviously leads to SU(2)l symmetry 
of weak interaction. The fermion field is a left handed doublet. Then, the 
current is of the form 



W 



e)L7^- 



(1.14) 



where r\ is the isospin triplet. Since the weak neutral current, J , has 



a right handed component, the neutral component of ( |1.14 ) is not a weak 
current. However, this can be identified with the electromagnetic current 



Jm = e h 



e^7„Q^ 



11-15) 



where Q is the charge operator with eigenvalue -1 for electron and the 
generator of £7(1) symmetry. Under orthogonal combination 

1 



h 



J 3 + - J Y 



(1.16) 



The consequence of the above equation is the relation, Q = T^ + 1V/2, 
weak counterpart of (|0|), where T^ is the weak isospin and Y\y the weak 
hypercharge. Like Q, Yyy generates the U(l)y symmetry group. Thus, this 
yields the SU(2)l x U(1)y gauge group. The fundamental fermions, quarks 
and leptons, come in three generations or families with weak isospin and 



weak hypercharge quantum numbers as shown in the Table 1.3. 



Table 1.3: Assignment of quantum numbers in SU(3)l X U(l)y 



generation 


n 


Yl 


1r 


Y R 


u c t 


1/2 


1/6 





2/3 


d s b 


-1/2 


1/6 





-1/3 


V e Vp v T 


1/2 


-1/2 








e \x r 


-1/2 


-1/2 





-1 



Electroweak Lagrangian 

The electroweak Lagrangian of the non-Abelian local gauge SU(2)l x U(1)y 
symmetry is of the form 

C ew = -\w i >^Wi v --B^B lxv + 2Y,¥ Ll ^¥ L (1.17) 
D^ = d^ + igfW^ + ig'YB^ (1.18) 

where g and g' are coupling constants. We get the photon field and the Z° 
field as the linear combination of W^ and B coupled to the electromagnetic 
current: 

Ay, = - sh\9 w Wl + cos 9 W B^, (1.19) 

Z^ = oxB w Wl + ^D.e w B ll (1.20) 

where tw.6w = g'/g, e = g sin 9\y The interaction part is 

Ci = - l-^JtW-^ + -^J~W + '^ + eJf A"l (1.21) 

l2V2 ** 2^/2 M M J 

The currents are given by 

J M = E^(l-75)VbifM9 + E^(l-75)« (1-22) 

^ c = E*-r>-«76)* ( L23 ) 

J e ™ = E*>Q* (1-24) 

where Vckm is the CKM matrix of flavour mixing (see below) and v = 
If - 2Q sin 2 W , a = Tf. 

Spontaneous breaking of symmetry 



In the electroweak Lagrangian in ( |1.18 ), introducing a mass term for both 



fermions and gauge bosons is forbidden by gauge invariance. However, the 
Higgs mechanism of spontaneous breaking of the SU(2)l x U(1)y symmetry 



[17| gives mass for the fermions and the three bosons of weak interaction. 



For the Higgs doublet field, $, 

* = ( to + ) (1-25) 



the Lagrangian is of the form 

C+ix) = (Mjtp,,*) - £/ x xi$x« + /» 2 ^ - A^*)' (1.26) 

x 

D^ = (l(d^ + i 9 -B^)+ig' T -±W^ (1.27) 

In ( |1.26| ), the three terms stand for free, Higgs-fermion interaction and po- 
tential parts as in the order given. Minimising the potential, ground state 
Higgs field is 

-(U)- -i/? (L28) 

This non-vanishing field describes the spontaneous breaking of the elec- 
troweak symmetry. Thus, the masses of fermions and bosons are generated: 

vg x 



C mass = - —Y^faaa+i^-j-j W+W 



•- / j 2 n J \ l w 3 ^ 



9 -99 

-gg' g 2 



+ -(W- B,) ( y _, J" ^ (1.29) 



where a stands for all the fermions, excluding the neutrinos. 

Masses for the bosons and the fermions (but not neutrinos) are, then, 
given by 

v 
m a = — f a , a = u,d, e, ... (1.30) 



M W = \g\ M z = y g z+g' 2 (1.31) 

M 7 = (1.32) 

Quarks mixing 

Mixing of different flavours is caused by the charged current mediated by 
W . The quark-mixing matrix, Vckm, hi ( 1-22 ) is a unitary (3 x 3) matrix 







/ Vud Vus Vub 



VcKM 



Vcd Vcs Vcb J (1-33) 

V v td v ts v tb 

By the above equation, the weak basis, q' and mass basis q are related by 
q'i = Vijqj The unitarity of the matrix V is: Vij — > exp(i(pi)Vijexp(—i(pj). 



9 



Further, the elements of the CKM matrix satisfy unitarity constraints: 
any pair of rows or columns are orthogonal. There are six constraints ( |1,33| ). 
For example, 

V ud V t * d + V us V t * s + V ub V t t = (1.34) 

This is useful to determine one element from the others. 

1.1.4 Weak Decays of Heavy Hadrons 

Knowledge of the weak decays of beauty hadrons, especially B mesons, is 
an important requirement to test the Standard Model fl9|| . Of about twenty 
input parameters, nine are the CKM matrix elements. The unitarity test 
of the CKM matrix demands precise prediction of its elements. The deter- 
mination of V u b, V c b, Vfb, Vt s , and Vtd, involves the b hadron weak processes: 
the charged current mediated modes, viz., the dominant b — > c + Y, and the 
CKM suppressed, b — > u + Y determine V c b and V u b respectively; the rare 
decays of the neutral current, b — > s + Y and b — > d + Y respectively de- 
termine Vts and V t d', and the weak decay, mediated by the charged current, 
t — > bW, determines V t b- (Y denotes either a lepton pair or a quark pair or a 
photon/gluon.). Precise determination of the elements of the CKM matrix 
is stringent need for checking the CKM model of the CP violation in the 
SM @. 

The heavy hadrons, hadrons containing a heavy quark whose mass itiq 3> 
A-qcd, is easier to study from the QCD point of view. For, the coupling 
constant is smaller at rriQ and the Compton wavelength of the heavy quark, 
1/rriQ, is much smaller than Aq CD . Thus, a lot of simplification arises in de- 
scribing the properties and decays of these hadrons. In the limit, tuq — > oo, 
with the resultant so called heavy quark symmetry, we get a low energy 
description of these heavy hadrons. This enables a model independent pre- 
diction of the CKM matrix elements \V u b\ and \V c b\. 
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1.2 Hadrons with Spin-flavour Symmetry 

In this section, we discuss the spin-flavour symmetry or heavy quark sym- 
metry exhibited by the hadrons containing a heavy quark, quark of mass 
rriQ 3> Aqcd- Then, we describe the theory of inclusive decays of heavy 
hadrons [g, [H], |2|, [23 1. 

1.2.1 Spin-flavour symmetry 

Because the heavy quark mass can be taken to the asymptotic limit, the 
study of weak interaction of the heavy hadrons becomes simplified 0, EI3, 



26, 27]: 



• the asymptotic freedom: the coupling constant is small enough to allow 
perturbative calculations; 

• the Compton wavelength of the heavy quark, 1/mn, is much greater 
than the typical hadron size, 1/Aqcd- 

In the limit tuq — » oo, the heavy quark acts as a static source of colour 
inside the hadron. Since the light degrees of freedom cannot resolve the 
distances smaller that the typical hadron size, the light stuff interaction with 
the heavy quark vanishes. Consequently, either state of the heavy quark spin 
is independent of the other. Thus heavy quark spin generates the SU{2) spin 
symmetry. On the other hand, the replacement of an infinitely heavy quark 
by another infinitely heavy quark does not change the description of the 
system. This leads to the flavour symmetry. Hence, for Nf heavy flavours, 
the hadron concerned is described by SU{2Nf) symmetry, called spin-flavour 
or heavy quark symmetry. 

In the static limit of the heavy quark, the four-velocity of the heavy 
quark is Vn = (1,0), i. e., the heavy quark is in the rest frame of the 
hadron containing it. Thus, both the decaying hadron and the corresponding 
final state are looking the same in this limit, since the product of the four- 
velocities of the initial and final state heavy quarks is unity. 

Also the third components of the heavy quark spin, ±1/2 and that of the 
light quark(s), j are separately conserved, J = j±l/2. With the heavy quark 
symmetry, the heavy hadrons can be classified by the quantum numbers of 
the light degrees of freedom, viz., flavour, spin, parity, etc., but independent 
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of heavy quark spin and mass pqj . As will be described in more detail later, 
the mass of the hadron is given by an expansion in the inverse powers of the 
heavy quark mass: 

M H = m Q + A + 0(l/m Q ) (1.35) 

where A stands for the energy of the light degrees of freedom. The term 
0{\/mo) is due to the kinetic energy of the heavy quark inside the hadron 
and the hyperfine mass splitting due to the spin orientation of the heavy 
quark with the light degrees of freedom. Its magnitude depends on the 
heavy quark mass, of the order of Aqqd- 

Under the breakdown of the SU(3) flavour symmetry, the c and b meson 
mass differences are of approximately same magnitude: 

m Ba - m Bd = A s - Ad (1.36) 

m Ds - m Dd = A s - A d (1.37) 

where the mass differences are respectively given by 90 ± 3 and 99 ± 1 MeV. 
Similarly, the mass difference between the baryon and the meson is given 
by: 

m\ b - m B = A Ab - A B (1-38) 

m Ac -m D = A Ac - A D (1.39) 

Their values are predicted to be 346±6 and 416± 1 MeV respectively. These 
mass differences receive corrections of 0(Aqc*£)/?71q). 

1.2.2 Effective Theory 

In the limit, ttiq — > oo, the QCD Lagrangian turns out to be a low energy 
effective one as the heavy particle is irrelevant at low energies [29, 30, [H], [3^, 



33, f3J, 35]. Now the theory is similar to the Fermi theory of weak interaction 
discussed in the previous section. 

As mentioned earlier, in the heavy quark limit, the heavy quark is in the 
rest frame of the parent hadron and it is almost on-shell. Its momentum is 

p^ = mQV » + k" (1.40) 



12 



where k is the residual momentum of O(Aqcd)- Accordingly, we define the 
propagator and the vertex function as 

' Y + f and ig s v a T a (1.41) 



Pq ~ m Q v - k 2 

respectively. In the above equation, (1 ± f)j2 are the positive and negative 
projection operators. They act finitely on the large and small components 
of the fields respectively: 

h v (x) = e im Q v - x P+Q{x), H v (x) = e- im Q v - x P-Q{x) (1.42) 

such that 

Q{x) = e- im Q v - x [h v {x) + H v [x)\ (1.43) 

Because of the projection operator, the field defined above satisfies v/h v = h v 
and \/)H v = —H v , where h v annihilates a quark of velocity v and H v creates 
an heavy antiquark. 

With these new fields, the QCD Lagrangian 

Cqcd = Q{x)(P - m Q )Q{x) (1.44) 

turns out to be of the following form for a heavy quark 

C Q = h v iv.Dh v - H v (iv.D + 2m Q )H v + h v (i P±)H V + H(i P±)h v (1.45) 

where Dj_ = D^ — v^, v.D satisfying v.D± = in the rest frame, the h v term 
corresponds to the massless degrees of freedom and the H v term to the heavy 
degrees of freedom which have to be integrated out, with fluctuation being 
twice the heavy quark mass. By using the equation of motion, the heavy 
degrees of freedom, H v , can be eliminated. That is by taking derivative of 
the above equation with respect to H v , we get the differential equation and 
its solution as 

(iv.D + 2m Q )H v = i p ± h v , H v = —-i p ± h v (1.46) 

IrriQ + iv. D 

Thus the nonlocal effective action results in 

Lett = h v iv.Dh v + h v i p±- — — —i P±h v (1-47) 

2m,Q + iv. D 
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where the second term corresponds to virtual processes. This Lagrangian 
can be expanded in powers of I/tuq: 

1 °° / ' J~)\ 

Ceff = h v iv.Dh v + V h v i p ± ^— i p ± h v (1.48) 

The Lagrangian is explicitly given as 
C eff = h v iv.Dh v + -^-h v (iD ± ) 2 h v + -^h v a flu G flu h v + 0(l/m 2 Q ) (1.49) 

Thus, in the limit txiq — > oo, only the first term survives, which is indepen- 
dendent of the heavy quark mass: 

£>hqet = h v iv.Dh v (1.50) 



For this Lagrangian, the Feynman rules are given already in (1.41|). In the 



above, we have applied the heavy quark limit to the effective Lagrangian, not 
the QCD Lagrangian. Thus, the heavy quark symmetry is not the symmetry 
of QCD but that of the low energy effective theory, (1.50|). The corrections 



appear to Chqet due to the finite mass of the heavy quark. As the heavy 
quark mass increases, the expansion becomes more and more powerful. 

In the rest frame, the 0(l/rriQ) operators are kinetic and chromomag- 
netic operators. 

Okin = -z h v (iD±) 2 h v -> h v (iD) 2 h v (1-51) 

2m Q 2itiq 

which is the gauge covariant extension of the kinetic energy due to heavy 
quark motion. 

Oman = ~r —h v (T uvG* 1 ' h v -> —h b SB c h v (1-52) 

4m Q r rriQ 

is the non-Abelian analogue of the Pauli interaction describing the colour 
magnetic coupling of the heavy quark spin to the gluon field. In the above 
equation for the magnetic operator, the spin operator S and the component 
of the colour magnetic field B l c are given by 

S i = -(^ ° t j, [S i ,S j ]=ie ijk S k (1.53) 

B l c = -- e ^ k G jk (1.54) 
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The corrections appearing at the next-to-leading order and beyond can 
be expressed as the short and long distance parts entangled. While the short 
distance physics is contained in the Wilson coefficients, the long distance 
physics are given by the local operators. 

As we have mentioned earlier ( 1.35| ), the energy of the light degrees of 



freedom A varies with the heavy quark mass. But the choice of the heavy 
quark mass is not unique. The mass previously referred to in this section is 
the mass in the Lagrangian. If one chooses it to be the pole mass, mo, then 
A is given in the asymptotic limit of the heavy quark mass as 

A ={M H -m Q ) | mQ _oo (1.55) 

If we define the heavy quark mass in some other scheme and at a certain 
scale, then A changes accordingly. However, if defined alternately as 

A - 5m = (M H - m Q ) (1.56) 

then the quantity A — 5m is invariant, where 5m is known as the residual 



mass of the order of Aqcd [36 1 . In that case, since the HQET allows the 



residual mass in the leading Lagrangian, Chqet would be of the form 

C-hqet = h v iv.Dh v - 5mh v h v (1-57) 

However, if the heavy quark mass is defined as the pole mass, then the 
residual mass term, 5m = 0. Now, these three quantities thq, A and 5m are 
non-perturbative parameters of HQET. They cannot be defined precisely in 
perturbation theory where ambiguities are of the order of Aqcd given by 
the renormalon graphs. 

1.2.3 Exclusive decays 

The heavy-quark symmetry nicely describes the exclusive weak decays by a 
single function |25|, |2fJ, |27j . With a mass- independent normalisation of the 
meson states, 

(M(p')\M(p)) = 2(p°/M M )(27Tf5 3 (p ~ P) (1-58) 

Thus, the meson state is redefined as 

\M(v) >= M^\M(p) > . (1.59) 
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Now, the meson state is characterised by the configuration of the light de- 
grees of freedom and this is the eigenstate of the effective Lagrangian which 
is supplemented by the QCD Lagrangian for the light degrees of freedom. 

Thus weak decays of a heavy meson state, in the heavy quark limit, 
corresponds to the replacement of the heavy quark, Q(v), in the initial state 
by another Q'(v'). The process is described by 

(M(i/)\h v 'fh v \P(v)) = £(v.i/)(v + v'Y (1-60) 

where the form factor £(v.v'), is independent of the heavy quark mass. 

Now, we see the weak decays of the heavy hadrons. There are six form 
factors associated with transitions of one flavour to another. The heavy 
quark symmetry relates them to one unique and universal function, known 



as Isgur-Wise function [27], £(v.v' = 1), where v,v' are the initial and final 
four- velocities. For the transitions due to the currents, Vn,, Au, we have 

(D*{v')\A, L \B{v)) = h Al (y)(y + l)e* 

-{h M v^ + h Aa {y)^V-v (1.61) 

(D(v')\V^\B(v)) = h+(y)(v+v'r + h-(y)(v-v'r (1-62) 

(D*{v')\V,\B{v)) = ih v (y)e^ a ^:v' a v p (1.63) 

The six form factors, in the heavy quark limit, are related to one another 
by £{y = v.v'): 

h+(y) = K(y) = h Al (w) = h A3 (w) = £(y), 

h-(y) = h A . 3 (y) = (1.64) 

Now consider the weak decay of B into D which amounts to the re- 
placement of the b quark in the initial hadron by the c quark. With the 
application of heavy quark symmetry, we have 

{D(v')\c vn »b v \B(v)) = t(y)(v + v'r (1-65) 



VM B M D 

The determination of the matrix elements by £(y) is facilitated by compar- 
ison with 

(D(v')\c v/ ^b v \B(v)) = f + (q 2 )(p+p'r ~ f-iP-p'T- (1-66) 
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Thus the independent form factors f± are given by £(y). 
M B ±M d 



f±(q' 



-ay), g 2 



M 2 B + M 2 D 



2yM B M D . 



;i.67) 



2y/ M B M D 

If £(y) = 1) the form factors are given by the hadron mass. Thus, the 
differential decay rate of B — ► -DZf is given by |37]] 

g^^ = ^3 IK*| 2 (M B + M D ) 2 M|>( y 2 - l) 3 / 2 ^) (1.68) 

The symmetry braking corrections, arising in the preasymptotic limit of 
the heavy quark mass, are 0(Aqcd/Mq) 2 - However, the matrix elements 
describing the 1/mg corrections vanish at the equal velocity limit. This is 
known as Luke's theorem J39[ |, an analogue of the Ademello-Gatto theorem 
[40|. The corrections of 0(1/toq) are studied in [38, 37, 41] 

Application of the formalism in which the weak exclusive decays are 
expressed in terms of the Isgur-Wise function is to extract the elements of the 



CKM matrix. From (1.68), we can determine \V c b\ in a model independent 



way m, g|, ||, ||, g§ |7], i 



51, §2, 53, 54, 7S]. 



1.2.4 Inclusive decays 

The theory of the inclusive decays of the heavy hadrons is formulated using 
the optical theorem [|j| |5g, |57], |58|, |59|, gO], q|, |6^, 63, Q. It allows one to 
determine the quantities like the CKM matrix elements cleanly in a model 
independent way. In this mode, considering the final state as the sum of 
all the hadronic states, we eliminate the bound state effects which are not 
known from the first principles of QCD. This is based on the idea of quark- 
hadron duality: the hadronic quantities can be calculated in terms of quarks 
and gluons at a particular kinematic point. That is, the hadronic quanti- 
ties specified in the physical region are calculable by means of a smearing 
mechanism in terms of quarks/gluons |3q] . 

Inclusive decay rate is expressed using the optical theorem as 

1 



T(H b - X) 



Mm 



-Im{H b \T\H b ) 



where T is the transition operator given by 



d 4 xT{C eff (x),C eff (0)} 



(1.69) 



;i.7o) 
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The effective Lagrangian C e ff is given by 

C eff = -2V2G f j£ c J CC , p. (1.71) 

Insertion of a complete set of states between the time ordered product, the 
decay rate is of the form 

F(H b - X) = _L£(27r) 4 <5 4 (pH - p x )\ (H b \C eff \H b ) | 2 (1.72) 

H x 



The Lagrangian ( 1.7l| ) receives a short distance correction due to the virtu- 



alities between Mw and m b . For the transition operator in ( [1.70 ), the OPE 



can be constructed in terms of a series of local operators containing heavy 
quark fields. These operators are arranged according to their dimensions. 

• The first operator is of dimensions three (D = 3): QQ (see Figure 
0a). 

• Secondly, D = 4 operator QiDQ is reduced to rriQQQ by the equation 
of motion (see Figur e|i~T| b). 

• Then the operator of D = 5 which is distinct from the D = 3 operator 
is the one involving the gluon field: Qgcr^G^Q 

Then the inclusive decay rate is given by these set of operators by 

(1.73) 
where c n are the short distance Wilson coefficients calculable in perturbation 
theory, and 

(0} H = -^ (H Q \0\H q ) (1.74) 

The hadronic matrix elements are expanded in inverse powers of the heavy 
quark mass (Figure 1.1) as 

(QQ) H = l- A{H) ~f G{H) +0{l/m%) (1.75) 
{Qg<r^CTQ) H = 2 l i%{H) + 0{l/m 2 Q ) (1.76) 
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Figure 1.1: Perturbative corrections to the Transition operator, T, (left) and 
to the local operators in OPE (right), (a) QQ (b) QgoGQ (Dotted lines are 
gluons). 
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where the HQET matrix elements are defined as 

AiH) = ^(H(v)\Q(zD) 2 Q\H(v)) (1.77) 

2/rrN 1 / TT,..\\*9 



Ag(H) = ^[H(v)\Q^ u G^Q\H(v)j (1.78) 

Hence, we have the decay rate of the form 

(1.79) 
In the above equation, the leading order term corresponds to the free 
quark decay rate. At 0(1/toq), the term /i 2 (H) describes the motion of the 
heavy quark inside the hadron and the term /j,q(H) the hyperfine interaction. 
The kinetic term for a B meson is estimated to be in the range of 0.1 to 0.5 
GeV 2 IH 67 , |68], fF?]] and the latest value is given by 



nl(B) = -Ai = 0.3 ± 0.2 GeV 2 (1.80) 

For the A& baryon, it is obtained using the mass values of the hadrons from 
the following relation 

(m Ab - m Ac ) - {m B - m D ) = \p*(B) - /4(A 6 )] (— - —J + 0(l/m 2 Q ) 

(1.81) 
where mB,D are the spin-averaged masses of B and D mesons. The kinetic 
energy of a heavy quark inside a meson and a baryon differs by 

nl{B) - AMb) = 0.01 ± 0.03 GeV 2 (1.82) 

It should be noted that the bound on the two terms are model independent 

fil(B) > n 2 G (B) w 0.35 GeV 2 (1.83) 

On the other hand, the hyperfine interaction term is determined pre- 
cisely from the difference between the masses of the vector and pseudoscalar 
mesons: 

lj, 2 G (B) = -(M B * - M B ) = 0.12 GeV 2 (1.84) 

For baryons it vanishes except in the case of the fig baryon. 



20 



Table 1.4: HQE predictions vs Experimental Values 





theory 


experiment 


t(S+)/t(B3) 

r(B° s )/r(B° d ) 
T(A b )/r(B°) 

Br(B) SL 
n c 


l+0(l/mg) 

(1.0 ± 0.01)+O(l/mf) 

0.98 + 0(l/mg) 

12.5% 

1.16 


1.02 ±0.04 

1.01 ±0.07 

0.78 ±0.05 

10.5% 

1.24 



The phase space factors corresponding to the C n are given by 



C3 = 


= 1 - 8x + 8x 3 - 


-x 4 - 


- Ylx 2 log x 2 


(1.85) 


C5 = 


- (I-*) 4 






(1.86) 



where x = m 2 /m 2 , the mass ratio of the final state quark to the decaying 



quark. The leading order factor given in (1.85) is exact [69] and the other 
due to higher orders in ( |1.86| ) is known only partially J73]. 

The total decay rate scales like m\. The leading order term corresponds 
to the free heavy quark decay rate. This is universal for all hadrons of a given 
flavour. The next-to-leading order terms appear at 0(l/m 2 ). They describe 
the heavy quark motion inside the hadron and the hyperfine interaction 
due to heavy quark spin projection. However, the next to leading order 
contributions are suppressed by m^ . 

With these, the important predictions for inclusive properties of beauty 



hadrons are given in Table 1.4. 



Because of the discrepancy in the case of A&, the semileptonic branching 
ratio of B mesons, Bt(B)sl |73|> 74, [T5| and the charm counting, n c , [72, [7q| , 
the third order in 1/mg terms have necessarily to be evaluated. This will 
be discussed in the next chapter. 
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Chapter 2 

Four-quark operators 



In this Chapter we discuss the evaluation of the expectation values of four- 
quark operators (EVFQO). The EVFQO are estimated in two different ap- 
proaches: (1) in terms of the form factor characterising the light scattering 
off the heavy quark inside the hadron, using a potential modef], and (2) 
from the difference in the total decay rates of beauty hadronsQ 

2.1 Preview 

As already mentioned, in order to systematically accommodate the contribu- 
tion coming from the terms at the third order of the heavy quark expansion, 
it is necessary to evaluate the expectation values of the four-quark operators. 
These are the dimension six operators involving both the heavy and light 
quark fields. Though their contributions are negligible as compared to the 
volume of heavy quark occupation, 0(Aqcd/'>tT'q) s , they are considerably 
enhanced due to partial compensation by the four-quark phase-space. As 
previously noted, the lifetimes of all mesons as well as baryons are almost 
equal among themselves respectively. But the knowledge of their individ- 
ual hierarchy depends on the outcome of the terms at 0(1/toq). Besides, 
the persisting confrontation of the theoretical estimate of the lifetime of 
Af, with the experimental value necessiates the evaluation of the four-quark 
operators. 

X S. Arunagiri, Colour-straight four-quark operators and lifetimes of b-flavoured hadrons 

(§• 

S. Arunagiri, A note on spectator effects and quark-hadron duality in inclusive beauty 
decays feol. 



22 



The hadronic matrix elements at 0{l/m%) is 

C(fx)(H b \(QT^Q)(qr^q)\H b ) (2.1) 

where T^ stands for some arbitrary Dirac structure, the coefficient C{n) 
describes the short distance QCD processes which will be discussed in the 
next chapter. The operators, in ([O]), are induced into a set of four by the 
spectator effects. They are: 



O q v _ A = Ml - 75)<797 M (1 - 7s)k 


(2.2) 


O q s _ P = 6(1-75)^(1 + 75)6 


(2.3) 


T$_ A = 67^1-75)^^(1-75)^ 


(2.4) 


T q s _ P = 6(1 - 75) A^(l + 75)^6 


(2.5) 


icuum states, we have 




(o^ A ) Bq = (ol- P ) Bq = fW B 


(2.6) 


( T v-*) Bq = ( Tq s-r) Bq =° 


(2.7) 



where fs is the leptonic decay constant of B meson. The evaluation of the 
operators is a difficult task. It involves the nonperturbative aspect of QCD. 
The expectation values of these operators are understood as the probability 
of finding both the quarks at the origin, ^(O)) 2 . The traditional way of 
estimating the expectation values of these operators for a meson is to use 
the vacuum insertion approximation and express the value in terms of the 
leptonic decay constant: 

^-(B\(QT,Q)(q^ q )\B) = |*(0)|| = *&£* (2.8) 

On the other hand, for baryon, one has to use the valence quark approxi- 
mation. These estimations are model dependent. We next discuss the eval- 
uation of the expectation values of four-quark operators made by Rosner 
PH , Neubert and Sachrajda |8^] and De Fazio and Colangelop^], followed 
by the evaluation made in potential model and from the difference in the 
decay rates of the SU(3) triplet mesons and baryons. 

In |8lJ , Rosner estimated the wavefunction density for A^ in constituent 
model using the DELPHI value of mass difference between E and S : 
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With the four-quark operators of a meson estimated in terms of the leptonic 
decay constant, the operators of a baryon can be made in a quark model 
using the hyperfine mass splitting of the 6-baryons. The mass splittings due 
to the hyperfine interaction in a meson, Mn, and in a baryon, Bijk, are 
given by 

AM(^) = ^a.^|*(0)|5 (2.9) 

J 9 rfiiVfij J 

AM(B tjk ) = ^a.£^3i>|*(0)|J (2.10) 

The value of (sQ.Sg) is (1/4, -3/4) for a ( 3 5i, So) Qq meson. With S qq = 1, 
(sQ.Sq) is (1/4, -1/2) for states with total spin (3/2, 1/2) for a Q(?g baryon. 
With these, the expectation values of four-quark operators for a 6-baryon is: 



For the baryon mass splitting, M(E£) - M(S 6 ) = 56 ± 16 MeV Q and 
for M(S*) - Af(S) = 46 MeV, we get 

|*(0)|i b = (2.6 ± 1.3) x 1(T 2 GeV 3 (2.12) 

As will be discussed in the next chapter, this value is found not sufficient 
to account for the smaller lifetime of A&. However, this work is in favour of 
larger value of |\I/(0)| 2 for Af, than that of B. 

Neubert and Sachrajda p2J parameterised the dimension six matrix el- 
ements in terms of hadronic parameters. The four-quark matrix elements 
between hadrons are given in terms of the hadronic parameters -Bj,ej. For 
B mesons 

' (B\(b L ^q L )(q L ^b L )\B) = ^-B 1 (2.13) 



2l7lB 



' (B\(b R q L )(q L b R )\B) = ?i^B 2 (2.14) 



2l7lB 



' (B\(b L ^t a q L ){q Ll n a b L )\B)^ f -^e l (2.15) 



2niB 



' (B\(b R t a q L )(q L t a b R )\B) = ^-e 2 (2.16) 



2rriB 
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where the hadronic parameters in the large- N c limit 

Bi = 0(1), et = 0(1/ N c ) (2.17) 

Using QCD sum rules, it is found that t\ = —0.15 and €2 = 0. 

Similarly, the matrix elements between A b states are expressed as 

' {^\{b Ll ,q L )(q L ^b L )\K b ) = -3^r (2.18) 



2m\ L 

where r = \^(0)\\ /\^(0)\%. It is expected that the ratio is about unity. 
But Rosner argued that it should be larger than unity. 

2.2 Potential Model Evaluation 

In the nonrelativistic quark theory, the wave function density and diquark 
density are related to the associated operator 

o 6 = (w*)(^rV), (2.19) 

where the r& )9 are arbitrary Dirac structures, through 
1 — - 2 



nj[g < B\(bb)(qq)\B > = \^>(0)\ 2 (2.20) 

JL- < A b \(bb)(qq)\A h > = J d 3 y\*(0,y)\ 2 (2.21) 

for a meson and a baryon respectively. The operators in ( |2.19| ) are colour 
singlet. The expectation values of these operators are related to the transi- 
tion amplitude of the light quark scattering off the heavy quark. Thus the 
determination is based on the knowledge of the light quark scattering form 
factor. 

The wave function at the origin, in momentum representation, is given 
by 

m = Jwr* {p) - (2 - 22) 

The transition amplitude is then the Fourier transform of the light quark 
density distribution: 

F(q) = ^- < H b (d)\qq(0)\H b (0) >= J d 3 x**(x)*(x)e^ x . (2.23) 
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Integrating over all q yields the expectation value: 

/ (^ F(q) = |X&(0)|2 = 2^ < H »\ lb ^(°)\ H b > (2-24) 

For any Dirac structure T, the light quark current density and the light 
quark transition amplitude are given by: 

J r (x) = qTq(x) (2.25) 

Mq) = ^ < H b(<l)\M0)\Hb(0) > (2.26) 

where the Jr(0) is gauge invariant operator, and it is not required to be a 
bilinear. Thus, for spin-singlet operators, we have 

< H(q)\bbJr(0)\H b (0) >= J ^_ < ^( q )| J r (0)|U 6 (0) > (2.27) 

And, similarly for spin-triplet operators, 

< H(q)\ba k bM0)\H b (0) >= J 7^ < S k H b (q)\J r (0)\H b (0) > (2.28) 
with S/2 being the b-quark spin operator and 

\S k H b >= f d 3 xba k b{x)\ H b > (2.29) 



Equations ( p.27 ) and ( p. 28 ) represent the general structure of four-quark 



operators. The above operators are local as required for by the HQE in the 
sense that the light quark operators enter at the same point as the heavy 
b-quark operators. These relations hold equally for different initial and final 
hadrons having different momenta smaller than m b . 

Equation ( |2.19 ) resolves into spin-singlet and spin-triplet operators for 



F b = 1 and T b = 775 (= cr) respectively. The light quark elastic scattering 
is described by the form factor F(q 2 ). In Ref. [85], the exponential ansatz 



and the two pole ansatz are used for the form factor. Both of them lead to 
a result which is the same for meson and baryon. 

Model for form factor 

As will be discussed in the following sections, the expectation values of the 
colour-straight operators are parameterised in terms of a single form factor 
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for both the B-meson and A& baryon [ p0[] . The extraction of the form factor 
involves an assumption of a function such that it satisfies the constraints on 
the form factor that F(q 2 = 0) is equal to the corresponding charge of the 
hadron. Then the form factor has to be extrapolated into the region where 
q 2 > 0. We take the hadronic wave function of ISGW harmonic oscillator 
model [ J36J I for the form factor. The wave functions of ISGW model are 
the eigenfunctions of the orbital angular momentum L = satisfying the 
overlapping integral which is chosen as the form factor. 

F(q 2 ) = N 2 exp[-q 2 /2(PJ + 2 )) (2.30) 

where N is the normalisation constant given by [2/3y /%/(/?? + (3f)] 3 ' 2 and /?'s 
are oscillator strengths. For the same initial and final hadrons (Pi = Pf = (3), 
the transition amplitude is [18 0| 



Ji V(g 2 ) = -?L = |*(0)| 2 (2.31) 



(2vr)3 ^ ' 4vr 3 / 2 

From the above equation, which is the master equation for evaluation of the 
EVFQO, it is obvious that the transition amplitude and hence the expecta- 
tion values of four fermion operators are proportional to the third power of 
the oscillator strength of the hadron. 

The calculation of /3's can be made using a QCD inspired potential. In 
the present calculation, we use the potential for a B-meson containing the 
Coulomb, confining and constant terms: 

V(r) = - + br + c (2.32) 

r 

where a = - 0.508, b = 0.182 GeV 2 and c = - 0.764 GeV. With quark masses 
m g = 0.3 GeV (treating m M = m^), m s = 0.5 GeV and m(, = 4.8 GeV, using 
a variational procedure with the wave function given in position space as, 

*(0)i s = ^e~^ 2 (2.33) 

the oscillator strengths of the beauty mesons are obtained. The values are 

(3 Bq = 0.4 GeV (2.34) 

l3 Bs = 0.44 GeV (2.35) 
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For the A& baryon, we use the same procedure as before but a potential 
of the form 

V(r) = -f^ + br + C'r 2 + c) (2.36) 

where £' = Cft 2 , with £ = 1 GeV and the r 2 term is a harmonic oscillator 
term justifying the consideration that the A& be a two body system using 
the same wave function as in ( |2.33 ), we determine the values of (3 for the 
baryons |8C]: 

(3 Ab = 0.72 GeV (2.37) 

/3 Sb = 0.76 GeV (2.38) 

We have taken the mass of the diquark system as 0.6 and 0.8 GeV for A& 
and E& respectively. 

The large value for the (3\ b is due to the presence of the 0(r 2 ) term in 
the potential. Otherwise, the value of /3a is no different than that of B s . 
These values are used in the subsequent calculations. A comment is in order 
on the choice of the same wave function for both the baryon as for meson: 
In the usual procedure, the ground state wave function for a baryon is 

* ground « e'^l^' 2 (2.39) 

where r\ tP are the internal coordinates for the three body system. Due to the 
idea of considering A& as a system containing the bound state of light quarks 
and a heavy quark, the separation between the two light quarks which make 
the bound state, is treated negligible. Then the baryon can be considered as 
a two body system, a reasonable approximation. The difference between a 
meson and a baryon is essentially due to the value of the oscillator strength. 

Expectation values 

We evaluate the expectation values of the colour-straight operators only for 
the vector and axial- vector currents. Nevertheless, the other currents can 
also be studied in the same fashion. Both the currents are possible for the 
B-meson, while axial currents vanish for A& baryon due to the light quarks 
constituting a spinless bound state. 

Essentially there is no difference between the exponential ansatz and the 
harmonic oscillator wave function in representing the behaviour of the form 
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factor, but they differ while fixing the scale: in the former case, the hadronic 
scale of one GeV is used whereas in the latter the same has been fixed by 
solving the Schrodinger equation. The two pole anstaz is based on the well 
founded experimental values. Basically, the use of the harmonic oscillator 
wave function of the constituent quark model is an alternate picture but in 
the very same lines of the ansatz. 

Hereinafter the operators are referred to by the following notation: for 
the meson 

< O q VA > H = (bT VtA bqr VtA q) (2.40) 

< T yA >H= (br V:A t a bqT V:A t a q) (2.41) 

where Ty, A = 7^,7^75 and t a t = 1/2 — 1/2N C . In what follows, q stands 
for u and d quarks and s for the strange quark. 

And, < Oy, A >, < Ty tA > will be respectively denoted as u>v, A , Ty, A - For 
the baryon, < 0(T)y > correspond to Ai,2 respectively. 

B-meson 

The parameterisation of the matrix element of the colour-straight operators 
for the vector current is 

< B(q)\J rv \B(0) >= -v^F B {q 2 ) (2.42) 

with the constraints F#(0) = 1 for valence quark current and F#(0) = for 
the sea quark current. The former is relevant for the b-meson composition 
of quarks bq. Then the corresponding transition amplitude is 

Under isospin SU(2) symmetry, 



A* v =< B(q)\O v \B(0) >=v ^J^F B (q 2 ) (2.43) 



For B s , we have 



(O v ) B - = 2.88 x 1(T 3 GeV 3 (2.44) 

(T V ) B _ = 9.61 x 1(T 4 GeV 3 (2.45) 



(O v ) Bs = 3.83 x 1(T 3 GeV 3 (2.46) 

{T v ) Bs = 1.28 x 1(T 3 GeV 3 (2.47) 
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If violation of isospin symmetry and SU(3) flavour symmetry is considered, 
then there comes Cabibbo mixing of eigenstates for d and s quarks among 
themselves. That is, d! = dcos9 c + ssin9 c and s' = scos6 c + dsin6 c . This 
would bring in negligible corrections. 

For the axial current, there are two form factors which are related to each 
other due to the conservation of the axial current, d^J^ = 0, in the chiral 



limit. By the Goldberger-Treiman relation [87] which equates axial charge 
form factor to the coupling gB*Bn at q 2 = 0, the operators involving the 
axial-currents are estimated in terms of a single form factor. Thus, given 
the value of the coupling g, the extraction of the transition amplitude is 
similar to the B-meson case. 

With the Goldberger-Treiman relation, we obtain the EVFQO for axial 
currents as 

G 1 (q 2 )=q 2 G (q 2 )=ge-« 2 / 4 ( 32 (2.48) 

where g = - 0.03 |p8[ p9|]. The numerical estimates are 



(0^> B _ = 8.63 x 10~ 5 GeV 3 (2.49) 

(T|) B _ = 2.88 x 10~ 5 GeV 3 (2.50) 

(O s A ) Bs = 1.15 x 10~ 4 GeV 3 (2.51) 

(T%) Bs = 3.84 x 10~ 5 GeV 3 (2.52) 

These values are insignificant as compared to the vector current contribu- 
tions. 

A5 baryon 

For the A& baryon, treating u and d quarks equally, 

{O q v ) Ab = 1.69 x 10~ 2 GeV 3 , ( T v) \ b = 5 - 64 x 10 ~ 3 Gey3 ( 2 - 53 ) 
In the case of E&, we have, 

( H °v) =2.01x10^ GeV 3 , / J2 T v) = 6.72xl0 -3 GeV 3 

\q'=u,d,s I a \q'=u,d,s I rj 

(2.54) 
Even though there are corrections to the form factors due to charge radius, 
they are ignored as we are interested in the wave function density at the 
origin. 
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Non-factorisable part of the FQO 

The nonfactorisable parts of the FQO come in four. The following is one of 
the ways of parameterising them 



— !p < B\(bq) V - A (qb) v „ A \B >= f B M B B l /2 (2.55) 

< B\(bt a q) v ^ A (qt a b) v _ A \B >= f B M B e x j2 (2.56) 



2M B 



< B\(bq) S -p(qb)s-p\B >= P B M B B 2 /2 (2.57) 



2Mb 
— ^ < B\(bt a q)s-p(Qt a b)s-p\B >= P B M B e 2 /2 (2.58) 

where B\ y2 and e\ )2 are hadronic parameters. They are related to Wv,A and 
Ty tA which are the expectation values of the operators Oy,A and Ty ;A as 
defined earlier: 

f B M B B 1 = 0! = A(t v + t a ) + — (w v + cja) (2.59) 

JVC 

/IMbSs = 2 = -2(ry - ta) - — (wy - w A ) (2.60) 

Die 



flM B e 1 = Pl = -^(rv + r A ) + (1 - ^ptXuv + wa) (2-61) 

1 , , 1, 1 



/iM B e 2 = p 2 = — (rv -r A )--(l- F )(wv - wx) (2.62) 



In the case of the A& baryon, the nonfactorisable piece corresponds to 

< (k)v-A(qb)v-A) = -^ A i - A 2 ( 2 - 63 ) 

2.3 Flavour dependence of decay rates 

In this section, we attempt to evaluate the expectation values of the four- 
quark operators from the differences in the total decay rates of the 6-hadrons. 
That means we assume that the heavy quark expansion, being asymptotic in 
nature, converges at 0(l/m 3 ) of the expansion. Already the next-to-leading 
order contribution due to terms of 0(l/m 2 ) is only about five percent of the 
leading order. Thus, it cannot be anticipated that the size of the terms at 
the third order in 1/m would be more than a few percent. On the other 
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hand, the observation made above is applicable only to the beauty case, 
since 

^C(„) (0 6 ) Hc » ^C(/z) (0,) Hb (2.64) 

" l c " l b 

where C(/j,) stands for some structure involving c± and {Oq) h , the dimension 
six FQO of hadron. Hence, in the background described, we make use of 
the total decay rates to obtain the expectation values of the four-quark 
operators for the 6-hadrons. Therefore, the present evaluation depends only 
on the heavy quark expansion and the SU(3) flavour symmetry, as has been 
shown in |Jl]] . 

Splitting of Decay Rates and Expectation Values of Four-quark 
Operators 

The B mesons, B~ , B° and B®, are triplet under SU(3)f flavour symmetry. 
The total decay rate of B meson splits up due to its light quark flavour 
dependence at the third order in the HQE. The differences in the decay 
rates of the triplet, T(B°) - r(fi"), T(B°) - T(B~) and F(B°) - T(B°), are 
related to the third order terms in 1/m by 

r(B )-r(B") = -r(,(l -xf {zi~(2c+ -c_ + 6) + 2(c + - c_/2 + 1)| 
x <0 6 ) B o_ B - (2.65) 

T(B°)-T(B-) = -r / (l-x) 2 |z 2 i(2 C+ -c_ + 6) + 2( C+ -c_/2 + l)| 
x (O 6 >fl0_ B - (2-66) 

r(B°)-T(B°) = -r' (l-x) 2 |(Z 1 -Z 2 )l(2c + -c_ + 6)} 

x (0 6 ) b o-b» (2-67) 

where r' = 2G^|y cfe | 2 m 2 /3vr, c_ = c^ 2 = {a(m)/a(m w )) 2 ^^ ,p = 11 - 
2n//3 (rif the number of flavours), 7 = 2, x = m^/m 2 and 

Zi = (cos 2 C (1 + |) + sin 2 9 c y/l - 4x(l - x) J (2.68) 

Z 2 = (sin 2 C (1 + |) + cos 2 <Wl - 4x(l - x)j (2.69) 



<Oe> = Q(bT,b)[(dT,d)-(uT, U )]j 



l(br il b)[(8T^)-(qr fl q)]) (2.70) 
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with q = u,d. In (2.65 - [037 ), the right hand side contains the terms cor- 



responding to the unsuppresed and suppressed nonleptonic decay rates and 
twice the semileptonic decay rates at the third order. 

For the decay rates T(B~) = 0.617 ps -1 , T(B°) = 0.637 ps" 1 and 
T(B°) = 0.645 ps" 1 @, the EVFQO are obtained for a B meson, as an 
average from ( p. 65 -2. 67): 



(O e ) B = 8.08 x 10~ 3 GeV 3 . (2.71) 

This is smaller than the one obtained in terms of the leptonic decay constant, 

IB- 

On the other hand, for the triplet baryons, A&, H~ and H°, with t(A&) 
< t(H°) ~ r(E — ), we have the relation between the difference in the total 
decay rates and the terms of 0(l/m 3 ) in the HQE, as 

r(A 6 ) - r(H°) = ^r' (- c+ (2e_ + c+ - 2) (o e ) Ab _ s0 (2.72) 

We obtain the EVFQO for the baryon 

(°6>a 6 -50 = 3.072 x 10~ 2 GeV 3 (2.73) 

where we have used the decay rates corresponding to the lifetimes 1.24 ps 
and 1.39 ps 0] of Aj, and H° respectively. The EVFQO for the baryon is 
about 3.8 times larger than that of the B meson. But it is about 3.5 - 4.0, 
depending on the changes in the mass and other sources of uncertainty. 

The expectation values for a baryon is quite large. In fact, Rosner [^] 
found it to be 1.8 times larger than that of a meson which accounts for 
about 45% (if hard gluon exchange is included) of the needed enhancement 
in the decay rate. However, what we obtain is model independent. We have 
only used the experimental total decay rates. The result is surprising as well 
as genuine, if one has to believe the experimental values used. Still, there 
is room for uncertainty of a few percent. On the other hand, the generic 
structure we employed can be decomposed into many as found in |53] by 
Neubert and Sachrajda. Using them, we will get an improved estimate, but 
in that case it may become somewhat model dependent. We should note 



here that preliminary study on lattice [92] shows that the EVFQO would 
enhance the decay rate of A&. 
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Chapter 3 

Spectator effects 



In this Chapter, we present the predictions for the ratio of lifetimes of beauty 
hadrons using the expectation values obtained in Ref. [80, |9(|. Also pre- 



sented is the determination the inclusive charmless semileptonic decay rate 
of At and discuss its influence on Br(b — ► X u lu) Q. 

3.1 Spectator Processes 



In the description of the inclusive decays of heavy hadrons by the heavy 

4- 



quark expansion, the total decay rate scales like m\. Corrections appear at 



0(l/m?) due to the heavy quark motion inside the hadron and the hyperfine 
interaction due to the heavy quark spin orientation. These corrections are 
about 10% of the leading order. As the heavy quark mass increases, the 
corrections are much suppressed by the inverse powers of the heavy quark 
mass. At the third order in l/rrc|, the effects are due to the participation 



of the light quarks (See Figure 3A). Though the effects are negligibly small, 



they distinguish the lifetimes of various hadrons of given flavour quantum 



number. We discuss first the processes involving light quarks [^, 94, ^, 
p6| . Then we will use the expectation values of the four-quark operators of 
mesons and baryons estimated in the previous chapter to get the lifetimes 
of the beauty hadrons. 



1 S. Arun agiri, Inclusive charmless semileptonic decay of At and Br(b — » X u lv), hep 



ph/0005266 
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Pauli interference (PI) 

For the decay of a (Qq) hadron, the quark level process is Q(qi) — ► Q'qfq'f(qi), 
where the subscript associated with g's denote the light quark in the initial 
hadron i and the light quarks produced in the decay /. If one of the final 
state light quarks happens to be the same as the initial state quark, then 
the wave functions will interfere. For example consider the decay mode 
B~ — > X c X q , where the subscripts c and q denote the charmed and light 
quark of either u and d final states. If either of the light quarks in the final 
state X q is identical to the one in the initial hadron, they will interfere de- 
structively. On the other hand, if Q' happens to be identical to %, then the 
interference will be a constructive one. The former will decrease the decay 
rate causing increase in the lifetime, while the latter will do the opposite. 

The constructive interference can happen only in the baryon. For exam- 
ple, in the charmless semileptonic processes of A^ — ► X u lv, the corresponding 
quark process is b(ud) — ► X u lv(ud). The u quark in the state X u and the 
u in the initial hadron will interfere constructively. This can happen in the 
nonleptonic process as well. 

VF-scattering (WS) and weak annihilation (WA) 

Another spectator process corresponds to the cross-channel of the quark 
decay. In the notation used before, the light quark in the initial hadron can 
scatter off the heavy quark through exchange of W field into a pair of light 
quarks, i. e., qiQ — ► q^q't- This process can happen only in baryons. This 
process is known as VF-scattering or IF-exchange, which occur in baryons 
only. 

Or, if the quark scattering off the heavy one happens to be an anti-quark, 
then they annihilate into a pair of leptons or quarks via a W boson: q^Q — ► 
qf q'f, where % is an anti-quark. This is the so-called weak annihilation 
process occuring only in mesons. 

3.2 Decay Rates due to Spectator Effects 

We summarise below the spectator contribution to the beauty hadron de- 
cays. 
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Figure 3.1: Spectator quark processes, (a) Pauli interference in B, (b) Weak 
annihilation in B, (c) W-scattering in A& and (d) Pauli interference in A&. 
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The quark level process b — ► cud: 

• B~(bu): the mZ interference term is given by 

Ar(fl-) = r 247r 2 Co< 0% > B - /ml (3.1) 



where 



C 4 - c 2 + (4 + ci)/7V c (3.2) 

7 



r = G 2 f m^|y XM | 2 /192vr 3 (3.3) 



c + =cl 1/2 = (a(m 6 )/a(m w )) 4 /* (3.4) 

with (3q = 11 — 2n//3,n/ being the number flavours. In (|3.1|) , (C^) 
stands for the EVFQO. 

• B®(bd): In this case, the (id annihilation contributes to the total decay 
rate. The magnitude of this effect is not significant. 

• B®(bs): No spectator effect occurs. 

The spectator effects become significant depending on the quark level 
process. 

For the above quark level process, in the case of beauty baryons, we 
have: 

• Ab(bud): There are two processes, dd destructive Pauli interference 
and uu W-scattering. The Pauli interference in this case enhances the 
decay rate resulting in the A& lifetime smaller than the average lifetime 



of the beauty hadrons. The combined decay rates are given in ( 3.14 
and 3.15J) in the latter part of this Chapter. 



• E°(bds),'E (bus): In the case of these hyperons, weak scattering and 
destructible Pauli interference occur. 

3.3 Lifetimes of Beauty Hadrons 

Using the results obtained for the expectation values of the four-quark op- 
erators of the hadrons, we present the results for the ratio of lifetimes of 
them. 
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The decay rates of the b-flavoured hadrons are given by 



T(H b - HM) = To 




(3.5) 


where, with x = m 2 /m 2 




f(x) = 1 - 8x + 8x 3 - x 4 + 12x 2 lnx 2 


(3.6) 


f'(x) = (1- 


-x) 4 


(3.7) 



are the QCD phase space factors and Ai and A2 correspond to the motion 
of the heavy quark inside the hadron and the chromomagnetic interaction 
respectively. We use: \ t (B) = 0.5 GeV 2 , Ai(A 6 ) = -0.43 GeV 2 and X 2 (B) = 
0.12 GeV 2 . Equation ( |3.5|) is further supplemented by the FQO at the order 
(1/m^) in the HQE. 

3.3.1 Ratio of Lifetimes of B and B^ 

Although the difference between the lifetimes of the charged and the neutral 
B-mesons is almost a settled issue, we check them once again using the expec- 
tation values of the colour-straight operators. This difference is attributed 
to PI and WA. Neglecting the WA as it is strongly CKM suppressed the 
result for the PI is 

ATf(B-) = ro247r 2 C <Q ^ >B '~ <Q ^ >B ' (3.8) 

where Co = c\ — c 2 _ + -^- (c+ + c 2 _ ) and the values for Wilson coefficients 

are: c+ = 0.84 and c = -1.42 with N c = 3. We use m b = 4.8 GeV and 

\V c b\ = 0.04. Then the ratio is 

This agrees well with the one obtained in terms of B-meson decay constants. 
The decay rates due to spectator quark(s) processes are: For B~, 

Ar»/(S-) = ( ^\V cb \ 2 (l - x) 2 [(2c 2 + - c 2 .)^ + 3(4 + c 2 _) Pl ] (3.10) 
Hence the ratio is 1.03. 
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3.3.2 Ratio of Lifetimes of B and B s 

The difference in lifetimes between the two neutral mesons B A , and B^ is due 
to W exchange. The numerical result is 

4|4 = 1-00 (3.11) 

Corresponding to the nonfactorisable part, we get the decay rate 

AT nf {B s ) = r (l - 4x) 1/2 i] (3.12) 

where Tq = G^m^/127r|V^b| 2 and rj is given by 

(2C+ " C ~ ] - ((1 " *Wi " (1 + 2s)#) + (C++ 2 C ' )2 ((1 - x)p{ - (1 + 2a 

(3.13) 



'/ 



3 |_ vv /ri v ,TZ ' 2 

The ratio becomes 1.02. 

3.3.3 Ratio of Lifetime of A b and B 

In the HQE, the difference in lifetimes between mesons and baryons begins 
to appear at order 1/m.Q. Nevertheless, it is dominant at third power in 
1/mQ. At this order, the FQO receives corrections due to WS and PI. They 
are 



T ws (A b ) = 92n 2 T cl < °^ 3 >Ab (3.14) 



r P/ (A 6 ) = -487r 2 r C 1 <Q ^3 >Ab (3.15) 

ml 

where C\ = — c+(2c_ — c+). As mentioned earlier, PI is destructible for 
radiative corrections and it enhances the decay rate leading to smaller life- 
time for Afc. The effect of WS, on the other hand, is colour enhanced and 
its consequence is smaller. Hence, 

^ = 0-79 (3.16) 

The decay rate modified by the nonfactorisable piece is given by 

Ar(A 6 ) = %^A[4(1 - x) 2 (c 2 - " 4) - t 1 " x ) 2 ( l + *0( c - ~ c +)( 5c + " c -)l 

(3.17) 
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Correspondingly, the ratio is 

r(A b ) 



0.84 (3.18) 



t(B-) 

In mesonic cases, the nonfactorisable piece gives a little higher value. In 
particular, the ratio of the lifetimes of the baryon and meson is significantly 
larger. 

3.4 Inclusive Charmless Semileptonic Decay 

Motivation 

Precise determination of the CKM matrix elements is fundamentally an 
important task for the Standard Model. There are as many as five of the 
nine elements of the CKM matrix that can be extracted from the knowledge 
of the weak decays of beauty hadrons. Particularly, \V cb \ and \V ub \ can be 
extracted in the framework of heavy quark expansion (HQE) in a direct and 
fairly model independent way |2l]] . With the measurements of the inclusive 
charmless semileptonic branching ratio of b hadrons by the ALEPH fl98[ , L3 



|H and DELPHI [|T0(| collaborations at LEP, the determination of \V ub \ 
has attracted renewed interest recently. The measured values of inclusive 
charmless semileptonic branching ratio of b hadrons and the corresponding 
values of \V u b\ extracted are 

= (1.73 ± 0.55 ± 0.55) x 1(T 3 , 

= (4.16 ± 1.02) x 10 -3 (ALEPH) (3.19) 

: (3.3±1.0±1.7)xl0~ 3 , 

: (6-o±?:oii:9±o.2)xio- 3 . (L3) (3.20) 

= (1.57 ± 0.35 ± 0.48 ± 0.27) x 10~ 3 , 

: (0.103iaoi2± °- 016 ±0.010). (DELPHI) (3.21) 



Br{b- 


-> XJvi) 




\v uh \ 


Br(b- 


-> XJui) 




\v ub \ 


Br(b- 


-> XJvi) 




\Vub/Vcb\ 



In these analyses, the measured quantity is 

(3.22) 



Br(b —* uev) V ub 2 



oc 



V cb 



Br(b — > ceu) 

This has an advantage of being free of many hadronic uncertainties that 
occur in the non-leptonic decays. Also, the determination of the CKM 
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matrix elements from inclusive decays can in general be made with less 
theoretical uncertainties than the ones extracted from exclusive modes. For 
semileptonic decays of baryons, however, there could still be large spectator 



effects due to Pauli interference as pointed out by Voloshin for H c [101], 
To the inclusive charmless semileptonic branching ratio of b hadrons, the 
contribution of A& is about 10% with the rest coming from the B mesons. 



There are many theoretical works [102] which study the B —> X u lvi to 



determine \V u b\. But not much has been done on the inclusive charmless 
semileptonic decay of A&. In the analyses of ALEPH and L3, A5 — ► X u lv\ 
is included in their sample, but in that of DELPHI, the rejection of kaons 
and protons that was used in the selection criteria reduces the contributions 
from A&. 

We study the inclusive charmless semileptonic decay of A^ |)7| . We find 
that the baryonic decay rate is larger by a factor of about 1.36, due to 
spectator effects, than that of the b quark decay rate. We then discuss the 
correction factors needed on Br{b — ► ulvi) to account for the spectator effect 
in A b ->• X u lvi. 

3.4.1 Inclusive Charmless Semileptonic Decay Rates 

The total rate of inclusively decaying beauty hadrons into a charmless final 
state is given by the HQEpT[ as 



G/IK,i.l 2 "» 5 / _ A t - 3A 2 _ 2 \2_ Q /_1 



where, with m being the mass of the heavy quark, the term of O(l/m ) cor- 
responds to the free heavy quark decay rate, the terms at 0(1/ m 2 ) describe 
the motion of the heavy quark inside the hadron (Ai = —0.5 GeV 2 for B and 
—0.43 GeV 2 for A&) and the chromomagnetic interaction due to the heavy 
quark spin projection (A2 = 0.12 GeV 2 ) which vanishes for baryons except 
Qq and the third order term in 1/m is given by 

C(fj,)(H\(br q )(qFb)\H). 

The operators in (...), denoted as (Oq) below, are evaluated in [|103[| by one 
of us for B and A&. The Wilson coefficients are describing the spectator 
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quark processes like Pauli interference, weak annihilation and VF-scattering 
which occurs in baryons only. 

In the mode A b — ► X u lvi, the spectator effect is the constructive interfer- 
ence of the u quark of the final state with the u quark in the initial hadron 
(see Fig. ??) which enhances considerably the decay rate. Otherwise, with- 
out spectator effects, the baryonic decay rate is almost the same as that of 
the heavy quark decay rate. The decay rate due to spectator processes is 
given by 

Ar(A b ) = 487r 2 r -^^ ( 3 - 24 ) 



where Tq = G 2 AV ub \ 2 7n 5 /l92-K 3 . Using the expectation values of the four- 
quark operators obtained in Ref. [103] for A b , (Oq) a = 6.75 x 10~ 2 GeV 3 , 



we obtain the ratio 

r(b -+ x u i n ) y ' 

whereas at 0(l/m 2 ), this ratio is 1.01. It should be pointed out that this en- 
hancement has a large theoretical uncertainty which is difficult to estimate. 
On the other hand, no such enhancement exists for A — ► X c eu e and thus 
the ratio is increased by 36% for A b . In passing, we note that the hadronic 
b — ► u decay of A& is not substantially enhanced due to the cancellations of 
constructive uu and destructive dd interferences. 

3.4.2 Br{b -> XJv) and V ub 

The inclusive charmless semileptonic branching ratio of b hadrons is given 
by 

Br(b - XJu t ) = (1 - f Ab )T(B -+ XJ^)t(B) + f Ab T(A b - X u l n ) T {k b ) 

(3.26) 
In principle, H & ' can also be produced at LEP. However, it is suppressed by 
ss production from vacuum relative to uu or dd production. Furthermore, 
the neutral Sj which consists of b, u, and s quarks has the same enhance- 
ment factor as A& in the charmless semileptonic decay rate, even though the 
charged H^ which is made of bds quarks does not have the corresponding 
enhancement. We will thus assume in this analysis that all weakly decaying 
b-baryons are A^'s. Thus using the above estimate with f\ b = 10%, t(B) 
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= 1.64 ps and t(A&) = 1.24 ps, we obtain the following value for inclusive 
charmless semileptonic branching ratio of b hadrons: 

Br(b -> XJui) = 1.16 x 1(T 3 (3.27) 

If the spectator effects are not included, then 

Br{b -f XJvi) = 1.13 x 10~ 3 (3.28) 



In the above estimate in Eqs. ( |3 . 2 5| - 3T28J ) , we have employed \V u b\ = 3.3 x 
10 -3 and m = 4.5 GeV. Also, we have used the decay rate for B mesons as 
estimated at 0(l/m 2 ) using the Eq. fl|!2| ), Br(B -> XJv h ) = 1.155 x 1CT 3 . 



No spectator effect occurs in b — ► ulv transition in B mesons except the 
negligible weak annihilation in B~ . Also, the experimental selection criteria 
used for the X u system mostly reject such final states. 

The corrections on \V u b\ due to the enhancement of A& — ► X u lv\ is about 
— 1.3%, and with these corrections, the central values for \V u b\ by ALEPH 
and L3 become, as given by \V u b\ = {Br expt / Brtheory) 0.0033, as 



4.03 x 10~ 3 (ALEPH) 
5.57 x 10~ 3 (L3) 



\ Vub \ ~ \ t;t;vv in-3 



Finally, we observe that (1) the contribution of A;, definitely influences 
the branching ratio due to a constructive Pauli interference and (2) such 
spectator effects in A^ decay have to be taken into account for a precise 
determination of \V u b\ which relies on charmless semileptonic decays. 
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Chapter 4 

Short distance 
nonperturbative corrections 



In this Chapter, we discuss the power corrections to the parton decay rate, 
due to renormalons, of the heavy hadron decays R. 

4.1 Renormalons 



As early as 1952, Dyson |104| showed that the perturbative expansion in a 
renormalisable theory like QED is not convergent. It is also true in QCD. 
This issue poses an important conceptual problem of giving a meaning to 
perturbation theory. Dyson's arguement is the following. An observable in 
QED, F(e 2 ) is expressed in an expansion in e 2 as 

F(e 2 ) = he 2 + / 2 e 4 + • • • (4.1) 

If there is a radius of convergence around e 2 = 0, then there will be a con- 
vergent result also for — e 2 . The latter corresponds to an attractive force of 
equal charges and a repulsive force of opposite charges. This would physi- 
cally mean that the state corresponding to the equal charges has unbounded 
negative energy. Thus the perturbative series cannot converge for negative 
e 2 . Hence the radius of convergence is zero. 

The divergent series can be dealt by the Borel summation technique 



1 S. Arunagiri, On the sho rt distance nonperturbative corrections in heavy quark ex- 



pansion, hep-ph/0009109 |103] 
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[ 105 , 106 ], A divergent perturbation series 

oo 

f[a] = Y / fna n = f 1 a + f 2 a 2 + --- 

n=l 

can be re-expressed by the Borel transform B[b\. 

oo ,„ ,2 

B[b] = £ /n+i-r « /i + /2& + /3 T + 
^o n! 2 

The original series F[a] can be reproduced by 



(4.2) 



f[a] 



db e~ b/a B[b] 



(4.3) 



(4.4) 



This would give a meaning to the series f[a] by the fact that the integral 
converges and B[b] has no singularities in the integration range. 
To be more explicit, let us consider 

d 4 x j** (0\T{j^x)j u (0)}\0) = (g^q 2 - q,q u )U(q 2 ) (4.5) 

With s = q 2 , the Adler function is defined as 

D(s) = -s-f-n(s) (4.6) 

as 

For the case of QED or QCD, we consider 

D n {Q 2 ) = J d A k F(k 2 ,kq,Q 2 )-^(blob) n (4.7) 

where we defined a = a(Q62)/-K and (blob) n represents n loops in a gluons 



line (Figure 4.1). Then the above equation would have the form 

k 2 ' " 



D n {Q 2 ) = J dk 2 k 2 <S> 



Q 2 k 2 



/3a log, Q2 



(4.8) 



where (3 is the first coefficient of the beta function. Finally, for momenta 
above and below Q 2 , we have the expressions corresponding respectively to 
the ultraviolet and infrared regions: 

Q 2 _ r 



D{Q 2 ) UV 



db e ~b/a_L_ 



D{Q 2 ) IR 



Q 2 

m — 1 



CO , 1 

b/a i_ 



db e 



B 



(4.9) 
(4.10) 



m— 1 
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Figure 4.1: Renormalon graph (dotted lines are gluons). 



where b = alog(Q 2 /k 2 )(h + 1) with h = m,n. 

In the infra red case, (|4.9| ), the poles in the Borel plane correspond 
to b = — (n + 1)//? on the negative axis. Hence the poles are not in the 
integration range. Thus the poles are at b = — 2/(3, — 3/(3, ■ ■ -. For the 
ultraviolet case, ( |4.10| ), the poles are at b = +1/(3, +2/(3, ■■ ■. In this case, 
b = (rn — l)alog(k 2 /Q 2 ). Thus, in this case, poles are in the range of 
integration and hence they are Borel summable. 

In the case of QCD which is asymptotically free in the UV region, the 
IR and UV renormalons defined above are interchanged. That is, the UV 
poles are at b = —1/(3, —2/(3, ■ ■ ■ and the IR poles at b = +2/(3, +3/(3, 

Operator Product Expansion 

The divergence of perturbation theory is deeply connected with the operator 
product expansion (OPE). Consider the OPE 

£>(Q 2 ) = [parton model] (1 + ai a(Q 2 ) + a 2 a(Q 2 ) 2 + ■■■)+ ^ C n (O n ) 

n=4,6,... 

(4.11) 
where C stands for the coefficient functions, (O n ) the operators of dimen- 
sion n. We cannot have a dimension two operators. So the term 0(1/Q 2 ) is 
absent. The perturbative part of the OPE receives the renormalon correc- 
tions |T2] . It has been shown that the absence of the first power suppressed 
term corresponds to the absence of singularity at b = +1/(3 in the Borel 
plane. However, the renormalon corrections to the perturbative series have 
one-to-one correspondence with the power suppressed terms. 
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4.2 Power Corrections 

The divergence of the perturbation theory at large order brings in an am- 
biguity to physical quantities specified at short distances. According to the 
present understanding, the ambiguity is given by a class of renormalon dia- 
grams which are chain of re-loops in a gluon line. The phenomenon is deeply 
connected with the operator product expansion (OPE). The perturbative 
part of the OPE receives the renormalon corrections fl2[| . Since in the OPE 
the first power-suppressed nonperturbative term is absent and the renor- 
malon corrections constitute short distance nonperturbative effect, they are 
more significant than large order corrections. 

The phenomenology of the power corrections is significant for the heavy 
quark expansion (HQE) to describe the inclusive decays of heavy hadrons 
by an expansion in the inverse powers of the heavy quark mass, tuq. As 
the inclusive decay rate of heavy hadrons scales like the fifth power of the 
heavy quark mass, the power corrections arise due to momenta smaller 
than the heavy quark mass. However, these IR renormalons would, be- 
ing nonperturbative effect, have greater influence in the HQE prediction of 
quantities of interest. These short distance nonperturbative effects can be 
sought for explaining the smaller lifetime of A&. We should note that these 
power corrections to heavy quark decay rate represents the breakdown of 
the quark-hadron duality. Therefore, it may shed light on the working of 
the assumption of quark-hadron duality in the heavy quark expansion. 

In this section, we study the renormalon corrections considering the 
heavy-light correlator in the QCD sum rules approach, assuming that the 
nonperturbative short distance corrections given by the gluon mass are much 
larger than the QCD scale. We carry out the analysis for both heavy meson 
and a heavy baryon. Our study shows that the short distance nonperturba- 
tive corrections to the baryon and the meson differ by a small amount which 
is significant for the smaller lifetime of the A&. 

Let us consider the correlator of hadronic currents J: 

n(Q 2 ) = i I (fixe iqx (0\T{J(x) J(0)}|0) (4.12) 

where Q 2 = — q 2 . The standard OPE is expressed as 

n(Q 2 ) « [ parton model] (1 + a x a + a 2 a 2 + ....) + 0(1/Q 4 ) (4.13) 
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where the power suppressed terms are quark and gluon operators. The 
perturbative series in the above equation can be rewritten as 

oo 

D(a) = 1 + a a + ^ a n a n (4.14) 

n=l 

where the term in the sum is considered to be the nonperturbative short 



distance quantity. It is studied by Chetyrkin et al [107| assuming that 
the short distance tachyonic gluon mass, A 2 , imitates the nonperturbative 
physics of the QCD. This, for the gluon propagator, means: 

D^ 2 ) = S J - <W (p + p) (4-15) 

The nonperturbative short distance corrections are argued to be the 1/Q 2 
correction in the OPE. 

Let us consider the assumption of the gluon mass A 2 3> A.q CD which is 
not necessarily to be tachyonic one. The feature of the assumption can be 
seen with the heavy quark potential 

Aa(r) 

V(r) = — + kr (4.16) 

where k ~ 0.2 GeV 2 , representing the string tension. It has been argued 
in |108| that the linear term can be replaced by a term of order r 2 . It is 
equivalent to replace k by a term describing the ultraviolet region. For the 
potential in ([4.16 ), 

k — ► constant x a\ (4-17) 

In replacing the coefficient of the term of 0(r) by A 2 , we make it consistent 



by the renormalisation factor. Thus the coefficient <r(A ) is given by [109]: 



18/11 

°W) = ^ k ' 2 ) [ 9^ ) ( 4 -!8) 




Introduction of A 2 brings in a small correction to the Coulombic term. By 
use of (4.18), we specify the effect at both the ultraviolet region and the 



region characterised by the QCD scale. Then, we rewrite ( f4.14| ) as 

D(a) = l + a a(l + ^ ] (4.19) 



T 2 
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where r is some scale relevant to the problem and k 2 should be read from 
(|4.18|) . We would apply this to the heavy light correlator in heavy quark 
effective theory. 

We should note that in the QCD sum rules approach, the scale involved 
in is given by the Borel variable which is about 0.5 GeV. But in the heavy 
quark expansion the relevant scale is the heavy quark mass, greater than the 
hadronic scale. Thus, there it turns out to be infrared renormalons effects. 
But, still it represents the short distance nonperturbative property, by virtue 
of the gluon mass being as high as the hadronic scale. 

Meson: For the heavy light current, J(x) = Q(x)ij^q(x), the QCD sum 
rules is already known flllOfl: 



fh 



-(A) = JL 



7T 2 Jo 



o^Ju?e- u l T D(a) 



+ 



1 



16r 2 



(gqaGq) + 



(4.20) 



where lu c is the duality interval, r the Borel variable and D(a) as defined 
in (4.14), but of the form defined in ( 4.1SJ ). It is, corresponding to the 
particular problem of heavy quarks, given as: 

-18/11" 



D{o)b = 1 + CLBOi 



1 A 2 

1 + — 



q(A 2 ) 
«(r 2 ) 



(4.21) 



where as = 17/3 + 47r 2 /9 — 41og(u;//i), with \i is chosen to be 1.3 GeV. 

With the duality interval of about 1.2-1.4 GeV which is little smaller 
than the onset of QCD which corresponds to 2 GeV and A > 0.6 GeV, we 

get 



A 2 = 0.35 GeV 2 . 



(4.22) 



Baryon: For the heavy baryon current 

j(x) = e abc (q 1 (x)C l5 tq 2 (x))Q(x) (4.23) 

where C is charge conjugate matrix, t the antisymmetric flavour matrix and 



a,b,c the colour indices, the QCD sum rules is given [112] by 



- f I e A/T 
2 lAb 



where 



1 



20vr 4 



6 



dLoco 5 e- w / T D(a) Ab +^E { 



n jo 



D(a) t 



o 
4tt 



OA h 



1 + 



7T 4 W 

(4.24) 
! (4.25) 
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with a Ab = rilog(2w/M) - r2). With /^ = 0.2 x 10~ 3 GeV 6 , (gg) = -0.24 3 
GeV 3 , (gqaGq) = m\ (qq), m 2 , = 0.8 GeV 2 , (aGG) = 0.04 GeV 4 and D(a) Ab 
is expressed in accordance with power correction factor found in [ |112| . As 
in the meson case, we obtain 

A 2 = 0.4 GeV 2 . (4.26) 

Now we turn to the heavy quark expansion. The total decay rate of a weakly 
decaying heavy hadron is, at the leading order, given by 



where 



1 - f (frfx) - i 



(4.27) 



r o - j -, nn _; - m (4.28) 



D IR = ~aa \ 1 + — ( -f^- ) | (4.29) 



Gj\V K M\ 2 ml 

1927T' 

As already mentioned, the power corrections are given by the IR renor 
malons: 

11/18N 

In ( 4.27] ), the factor £ corresponds to the IR renormalons which corresponds 
to the square root of the A 2 term in the above equation. These corrections 
are estimated to be O.iro and O.liro for B and A^ respectively. This is 
significant in view of the discrepancy between the lifetimes of B and A& 
being 0.2 ps" 1 with T(B) = 0.68 ps" 1 and r(A 6 ) = 0.85 ps" 1 . 
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Chapter 5 

Quark-hadron duality 



In this Chapter, we discuss qualitatively the isssue of the quark-hadron 
duality, for breveity duality. It is based on the assumption of the convergence 
of the heavy quark expansion (Chapter 3) and the corrections arising due 
to renormalons contributions (Chapter 4). 

Prologue 

The issue of duality is as old as QCD itself. The idea of duality can be 
formulated in many ways. The central idea stems from the fact that the 
hadronic quantities are defined in the Minkowsky region, whereas in terms 
of quarks and gluons they are specified in the Euclidean region. The two 
regions are connected by analytical continuation which is done by hand. 
However, the lore is that the hadronic observables can be calculated in terms 
of partons at a particular kinematic region. In their classic paper, Poggio, 
Quinn and Weinberg proposed a procedure, known as smearing, that the 
calculation of physical quantities can be obtained in terms of partons by 



averaging over a suitable energy range 65]. 

In the OPE procedure, as outlined in the previous Chapter, the hadronic 

quantities are given by the leading parton quantity plus some local operators 

representing nonperturbative aspects. For example, for the leptonic decay 

constant of a meson, say B, we have 
/-wo 
f B ~ / did p(uj) + nonperturbative corrections (5-1) 

Jo 

where p(u>) represents spectral density of quark(s) and u> is the duality inter- 
val. The determination of the duality interval depends on the contribution 
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coming from the nonperturbative physics. On the other hand, the pertur- 
bative corrections to the leading term is also important. However, for any 
quark level quantity, the perturbative corrections are not known beyond the 
first few orders. 

The inclusive decays of the heavy hadrons are described by the OPE 
based expansion of the hadronic matrix elements in the inverse powers of 
the heavy quark mass. The idea of quark-hadron duality is the underlying 
assumption. The discrepancy between theory and experiment for quantities 



given in Table 1.4 casts doubt on the validity of the assumption of duality. 
According to the HQE, the total decay rate of a heavy hadron scales like 
the fifth power of the heavy quark mass concerned. As we have seen in the 
previous Chapters, the corrections to the leading order (partonic) decay rate 
arises due to bound state effects which start appearing at the second order 
in 1/rriQ. These corrections are of the order of about 10% of the free heavy 
quark decay rate. It is significant for the HQE predictions. 

In this context, [ pL13| j have shown that the violation of duality in HQE is 
of exponential/oscillating in nature: 



n(Q 2 Uoiation = e~ CQ /A Qon (5.2) 

where C is constant and Q 2 is the energy scale. However, this violating 
effect is not quantified. They attribute this quantity for the discrepancy in 
the inclusive properties shown in Table (|1.4|). 



On the other hand, in [114], the weak decay of heavy hadrons is studied 
in the 't Hooft model. It has been found that the duality holds good with 
the presence of terms of order 1/rriQ. Such a term is absent in the HQE. We 
should note that the first-power-suppressed term is absent in the OPE itself. 
We should note that the 't Hooft QCD is 1+1 dimensional where confinement 
is bulit-in. But, in QCD, the phenomenon of confinement is not understood. 
Therefore, we cannot expect every aspect of the two-dimensional QCD to 
agree in toto with the QCD. 

With these remarks, we discuss briefly below the qualitative aspects 



derived from the study in Refs. [90, 103 1. 
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Perturbative corrections at large order 

The renormalon corrections evaluated in the previous Chapter for the B 
meson and A& baryons are obtained in terms of the gluon mass, A 2 = 0.35 



and 0.4 GeV respectively [103]. The assumption on A is made to repre- 



sent the short distance nonperturbative effect. They correspond to the IR 
renormalons in the description of the inclusive decays: 

-9/11 



AT(H Q ) IR « a a s x 



A 2 / a(X 2 



Numerically, the IR renormalon corrections are found to be 



(5.3) 



AT(B) IR « 0.ir (5.4) 

AT(A b ) IR w 0.1ir (5.5) 

where Tq 6-quark decay rate at the tree level. Since the sign of these contri- 
butions is fixed to be positive, the short distance nonperturbative corrections 
yield a small but significant enhancement of the A& decay rate. 

These corrections have to be construed as duality violating effects. Quan- 
titatively, the violations amount to be about 10%. In recent literature, these 
duality violating effects are argued to be the 1/Q 2 corrections. However, it 
cannot have operators representations. In the case of HQE, if it is considered 
to be 1/mn corrections, we cannot have the corresponding operators. 

Our conclusion is that these corrections are important in predicting the 
inclusive properties in HQE. 

HQE for beauty mass 

In the Chapter 2, we have evaluated the EVFQO from the differences in 
the inclusive decay rates [90]. Assuming the convergence of the heavy quark 



expansion is valid one as long as there is an uncertainty of few percent. 
This would also positively suggest that the quark-hadron duality violating 
oscillating component might be small. So such violation would not deter a 
decent determination of quantities of interest in the heavy quark expansion 
like the lifetimes of beauty hadrons, the semileptonic branching ratio and the 
charm counting. For the reasons mentioned in the beginning, the assumption 
on convergence cannot be made for charmed case. 
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There are renormalon contribution from the perturbative part of the ex- 
pansion. They are IR renormalons of the order Aqcd- We, as of now, don't 
have deep insight of it. They would be expected to differ for meson and 
baryon. Because this contribution does not correspond to any local opera- 
tors of the theory and independent of the heavy quark mass, a difference of 
about 50 to 100 MeV between meson and baryon would imply much signifi- 
cance for the quantities described by the heavy quark expansion. We cannot 
on obvious terms argue that renormalons are related to the assumption of 
quark-hadron duality. On the other hand, it will shed light on the quantities 
concerned and the underlying assumption versus the heavy quark expansion. 

Our studies show that the inclusive decays of heavy hadrons can be 
studied within the framework of the HQE, notwithstanding the aspects like 
exponential violation [115, 116] which have not been quantified. 
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Chapter 6 

Conclusion 



In this Thesis, we have addressed two important issues, namely, the evalua- 
tion of EVFQO and the quark-hadron duality in the heavy quark expansion. 
Our work primarily shows that the existing discrepancy between theory and 
experiment over the ratio of lifetimes of A& baryon and B meson can be ex- 
plained within the heavy quark expansion. The assumption of quark-hadron 
duality is found to be reasonably working well in the beauty hadrons, since 
the b quark is sufficiently heavy. 

The results of our study concerning the EVFQO and the ratio of lifetimes 
of Aft and B are given below: 

• The ratio of EVFQO of A& and B should be greater than unity, contrary 
to common expectations. Precisely, it should be in the range of three 
to four. Our estimation in the potential model is supported by out 
another evaluation using the difference in decay rates of SU(3) triplet 
hadrons. The latter is r model independent result, which depends only 
on the heavy quark expansion and the SU(3) flavour symmetry. Our 
result for the ratio t(A;,)/t(-B), from the potential model, is: 

r(A 6 )/r(B) « 0.79 - 0.84 (6.1) 

where the lower limit is fixed by the factorisable parts of the FQO and 
the upper limit by the inclusion of non-factorisable parts. The same 
quantity, in the model independent approach, is in the range of about 
0.78 to 0.81. 
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• The ratio |*(0)|^/|*(0)|g is about 3.8 from the potential model es- 
timation, and about 3.5 from the model independent approach. The 
latter value and the corresponding prediction would change, if the var- 
ious structures of FQO are included. In both cases, we have estimated 
the V — A structure only. For the purpose calculating the ratio of 
lifetimes, what we have done is sufficient. To calculate the absolute 
lifetime, one has to take into account all the possible structures of 
currents. 

• The corrections due to renormalons causes an increase in the decay 
rate of At,. This should also have to be taken into account in the 
calculation of the inclusive properties of the heavy hadrons. 

Our potential model, approach has an apparent draw back that the potential 
has a term of r 2 term for baryon. In fact, the potential that generates 
the heavy hadron mass spectrum is different. However, this potential too 
satisfies the Schrodinger equation. We have used the same wave function 
for meson and baryon. For the latter case, it is an approximation. 

Recently, in [B2], an exploratory study of the evaluation of the EVFQO 
of A& showed that the decay rate of the baryon can be significantly enhanced 
due to the third order term in the heavy quark expansion. On the other 
hand, in the QCD sum rules approach, the estimated EVFQO yields the 
ratio of lifetimes of A& and B to be 0.79. 

For the first time, we have studied the inclusive charmless semileptonic 
decay of A&. It is significant, in view of the experimental estimate of the 
inclusive charmless semileptonic branching ratio of beauty hadrons. This 
quantity is useful for a clean extraction the CKM matrix element, \V u b\. We 
have observed that in measuring Br(b — ► X u li>), one should have to take 
the spectators contribution in the baryon. Our prediction of branching ratio 
including spectator effects is: 

Br(b -> X u lv) = 1.16 x 1(T 3 (6.2) 

This result is significant in view of the recent measurement of these quanti- 
ties at LEP by ALEPH, L3 and DELPHI collaborations. 

Our observations on the issue of quark-hadron duality in the heavy quark 
expansion are the following: 
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• For the 6-quark mass being asymptotically heavy enough, the inclusive 
decay rate of beauty hadrons apparently converges at the third order 
in 1/777,5. But, it is not so in the case of the charm sector. As the results 
show, the EVFQO obtained for the B meson and A& baryon yields the 
ratio of their lifetimes closer to the experimental value. There is no 
reason to expect contribution to come from beyond the 0(l/m?). 

• The implication of the assumption of the convergence of the HQE is 
that quark-hadron duality works well. On the other hand, the expo- 
nential/oscillating violation of duality cannot be expected to be more 
than a few percent. 

• Our study of the renormalons contribution to the perturbative decay 
rate is important. If it is a duality breaking correction, as we esti- 
mated, it is about 20% of the free heavy quark decay rate. It may be 
taken as the 1/rrib correction, but we don't have an operator of appro- 
priate dimension. This means that we rule out the possibility of the 
absence of the I/ttt-q term being responsible for the duality violation. 

Thus our conclusion is that the quark-hadron duality holds good in the 
HQE. The violation of few percent is within the control. 

In view of the important task of predicting the elements of the CKM 
matrix, the clear idea of inclusive decays of heavy hadrons would be useful 
to reduce the error. Our studies are simple demonstration concerning the 
two aspects of the heavy quark expansion, namely, the systematic accom- 
modation of the spectator effects and the knowledge of the extent to which 
the notion of duality is reliable. 
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